Electron energy loss and induced photon emission in photonic crystals 
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The interaction of a fast electron with a photonic crystal is investigated by solving the Maxwell 
equations exactly for the external field provided by the electron in the presence of the crystal. The 
energy loss is obtained from the retarding force exerted on the electron by the induced electric field. 
The features of the energy loss spectra are shown to be related to the photonic band structure of the 
crystal. Two different regimes are discussed: for small lattice constants a relative to the wavelength 
of the associated electron excitations A, an effective medium theory can be used to describe the 
material; however, for a ~ A the photonic band structure plays an important role. Special attention 
is paid to the frequency gap regions in the latter case. 
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I. INTRODUCTION 

Near-field spectroscopy can be performed using 
spatially-resolved electron energy loss spectroscopy in 
scanning transmission electron microscopes, where the 
electron beam provides an external evanescent field with 
which to probe the sample within a spatial range on the 
nanometer scale. In particular, the relatively intense, 
low-energy part (< 50 eV) of the loss spectrum can be 
used with minimum sample damage to provide some in- 
sight on plasmons and other collective excitations. 1-4 
For the electron velocities typically employed in micro- 
scopes (above half the speed of light) and for samples that 
are homogeneous across distances of a few nanometers, 
frequency-dependent dielectric functions are sufficient to 
describe the materials that are involved and the loss spec- 
tra reflect the geometry of the sample interfaces. 5-10 

The geometry becomes particularly important in pho- 
tonic crystals, where the periodic spatial modulation of 
the dielectric function affects the propagation of radia- 
tion to the point of forbidding it within band-gap energy 
regions. As a consequence, photonic crystals are known 
to inhibit the spontaneous emission of light within the 
band gap. 11:12 They can also be used to make omnidi- 
rectional dielectric mirrors that reflect light from all di- 
rections without absorption 13 and wave guides able to 
deflect light around sharp corners on the scale of the 
wavelength. 14 These are applications of photonic crystals 
that involve free external radiation, but equally remark- 
able effects arc expected to accompany evanescent fields 
like those of external electrons. 

In this work, we examine the energy loss spectra of 
electrons moving near or inside photonic crystals over 
a wide range of lattice parameters a. When a is much 
smaller than the wavelength of the radiation associated 
to a given energy loss, A, it is shown that the crystal can 
be described to some extend by an effective dielectric 
function. 15-18 For larger lattice constants, the photonic 
band structure becomes more complicated and this is re- 
flected in the loss spectra. The relation between the band 
structure, the reflectance of photonic crystal slabs, and 



the electron energy loss spectra is discussed in detail in 
what follows. 

In a previous development, Pendry and Martin- 
Moreno 19 calculated the energy loss probability for elec- 
trons moving near crystals made of either metallic 
spheres or metallic cylinders in the a << A limit us- 
ing the transfer matrix approach 20 to solve Maxwell's 
equations. Their results exhibit a complicated energy- 
loss structure even for relatively dilute crystals , which 
has not been reproduced in the present work. Therefore, 
in order to test the convergence of our method, we have 
calculated the energy loss near the surface of a crystal 
of dilute spheres and have found very good agreement 
between the results derived from the theory described in 
this work and those obtained from the analytical expres- 
sion for isolated spheres. 9 ' 10 

In addition to producing energy loss, the interaction 
between the electron and the crystal gives rise to the 
emission of the so-called Smith-Purcell radiation. 21 This 
effect has already been discussed for one-dimensional 10 
and three-dimensional 22,23 crystals. Radiation emission 
is one of the contributions to the total energy loss, and 
in frequency regions where the crystal is transparent the 
probability of these two must coincide. Examples of this 
are offered below. 

Here, the energy loss and photon emission probability 
are calculated from the induced electric field, which is 
derived using the reflection coefficients of the crystal, 19 as 
explained in Sec. II. Results for the case of small lattice 
constants (a << A) are given in Sec. Ill and for larger 
ones (a ~ A) in Sees. IV and V. The main conclusions 
are summarized in Sec. VI. Gaussian atomic units (a.u., 
that is,?i = 77i = e = l) will be used from now on, unless 
otherwise specified. 



1 



II. THEORY 
A. The reflection coefficients of the crystal 

We shall consider crystals that are composed of a num- 
ber of layers perpendicular to the z direction and ex- 
tending from z = towards z < 0. The host region 
outside the crystal will be assumed to be described by 
a frequency-dependent dielectric function Ch{u>) and a 
magnetic permeability /Lt/j(w). Each crystal layer con- 
sists of the repetition of a given object with certain two- 
dimensional translational symmetry that is shared by all 
layers. The crystal will be characterized here by its reflec- 
tion coefficients, that is, the amplitudes of the reflected 
plane wave components for a given incident wave. 

In general, an external plane wave of frequency u> that 
propagates near the crystal can be represented as 

E± exp(iK± • r), 

where Q = (Q x , Q v ), 

K± = (Q,±ir Q ), 

r g = Q 2 ~ k h , 

and kh = (w/c)(e^ / u/ l ) 1 / 2 is the momentum of the plane 
wave. The real part of Tq is chosen to be positive, and 
the + (— ) sign in these expressions stands for a wave 
moving towards z > (z < 0). 24 

The momentum Q parallel to the surface has been sin- 
gled out to make explicit use of crystal symmetry: in- 
voking momentum conservation, a wave incident from 
the z > region with momentum Kq will only produce 

a discrete set of reflected waves of momentum Kq +G , 
where G runs over reciprocal surface lattice vectors. 

The transversal character of these waves (i.e., the fact 
that Kg - Eg = 0) can be exploited to express the electric 
field in terms of s and p components according to 

e q = E Q,s Cq,s + e q, p eq,p, 
where the vectors 

^Q,s = q(-Qv>Qx,0) 

and 

satisfy the identities s ■ eq g = eq • eq = 1 and 

eq, s ' eq iP - e± s • K± = e± p ■ K± = 0. Also, s is 
perpendicular to the plane denned by Q and the surface 
normal. 

The amplitudes of the reflected waves depend linearly 
on the amplitudes of s and p components of the incident 



wave, and the coefficients of the linear relation between 
them are the reflection coefficients -Rqo, implicitly de- 
fined by 

H +G j r = J2 r qg[Eq,*>y> (i) 

a' 

where a and a' run over polarization directions s and 
p, and the super-indices r and i stand for reflected and 
incident components, respectively. 

We have used the layer KKR method to calculate the 
reflection coefficients Rqq, both for a << A and for 
a ~ A. In the layer KKR method, the transmission and 
reflection coefficients are calculated exactly for each sin- 
gle layer using multiple scattering in a basis set of multi- 
poles centered around each object of the layer. Scattering 
among layers is then expressed in terms of those coeffi- 
cients. The maximum multipole order Z max and the num- 
ber of reflected and transmitted beams (i.e., the number 
of G vectors) are the only convergence parameters, which 
have been tested in all calculated results shown below. 
Stefanou et al 25 developed this method for spheres and 
we have extended it to be able to deal with non-spherical 
objects and arbitrary values of Z max . Further details of 
the method will be given elsewhere. 26 

B. The field of the electron 

We shall consider an electron moving with constant 
velocity v along a trajectory parallel to the crystal surface 
and described by r t = (vt, y$, z ), with z > 0, so that the 
electron moves in the host medium described by eh and 
/i/j. Neglecting the crystal for the moment, the electron 
electric field can be expressed in frequency space u> in 
terms of vector and scalar potentials as 

iu) 

E = — A - V0o- 

c 

In the Lorentz gauge, Maxwell's equations can be recast 
as 

(V 2 + fc 2 )^ = --P 
eh 

and 

(V 2 + fc 2 )A = -^j, 

where p(r, u) is the electron charge density, and j(r, u) = 
(vp, 0, 0) is its charge current. 
Using the relation 

<*+*>/g iid^ --*«M 

and noticing that p(r, t) — — S[r— (vt, yo, zo)], the electric 
field is found to be 



2 



Eo(r,a/) = [— 



rfp r 

27T2 J 



dte 



iutt 



e ip-[r-(vt,y ,z )] 

P 2 -k 2 h - iO+ ' 



where the time integral represents the inverse Fourier 
transform that permits obtaining p(r,oj) in terms of 
p(r, t) . The above integral can be reduced to 

Eo(r,o;)= /"d^e^-^^E* (2) 

where Q = (u>/v,Q y ) and 

w T Q ve h c 2 

The + (— ) sign must be used in these expressions when 
z > z (z < z ), so that the integrand of Eq. (2) rep- 
resents a plane wave that propagates towards positive 
(negative) z's. When the electron is moving in vacuum, 
Tq is real and the waves in the integrand of Eq. (2) are 
evanescent. 24 However, when kh is real and larger than 
lu/v, some of those waves describe Cherenkov radiation 
that propagates without attenuation; this is the case of 
electrons that travel faster than light in the medium. 



C. The field induced by interaction of the electron 
and the crystal 

The decomposition of these plane waves into s and p 
components is readily found to be 

[ ^ sl ~ qt q c 2 

and 

\ F ± n _ , fcfe -iK^-(0,yo, Zo ) 

[ ^ ~ ± vQe h ° 

Each of the incident plane waves [i.e., each value of Q y 
in Eq. (2)] gives rise to a set of reflected waves whose 
amplitudes are obtained from Eq. (1). Therefore, the 
electric held in the region near the ion can be constructed 
as the sum of E and the reflected held, 

E(r,w) = E (r,w) + E r (r,w), (3) 

where 

E r (r,w) = 

J2 I dQ v e iK 5 + c r R% G [E^J e+ +G ff , (4) 

and the integral over Q y has been copied directly from 
Eq. (2) in virtue of the linearity of Maxwell's equations. 



D. Electron energy loss in front of the crystal 

The electron energy loss can be calculated from the 
retarded force exerted by the induced part of the electric 
held E lnd acting back on the electron. Integrating this 
force along the trajectory and dividing by the total path 
length L, one finds 



AE 
~Ax 



jjdt 4 nd (r t , t)= f udu P(lu) 



where 



P(u>) 



TTLijL 



dtRe{e~ llJJt E^ d (r t ,uj)} 



(5) 



is the loss probability per unit of path length. 

For an electron moving parallel to a crystal surface, 
Eq. (3) permits separating the loss probability as 



P = P + P r , 



where 9 



I v 2 I a 2 - k 1 

PoM = —im{(-p h - -)io g [ - g ;. 2 h 



]} 



corresponds to the contribution of Eo in the absence of 
the crystal. Here, q c is a momentum cut-off related to en- 
ergy conservation. This contribution vanishes in vacuum. 
However, it gives rise to Cherenkov losses when eh and ph 
are real, in which case 27 ' 9 Po = |^h|[l/c 2 — l/(v 2 £hPh)], 
subject to the Cherenkov condition v 2 e h ph > c 2 . 

The remaining contribution P r is due to the retarding 
force exerted by the reflected field E r . Inserting Eq. (4) 
into Eq. (5), the time t and the lateral impact parameter 
yo appear only through the factor exp{i(G x vt + G y yo)}, 
so that time integration eliminates all vectors G with 
G x ^ from the sum in Eq. (4). Furthermore, averaging 
over yo leaves only the G y — term and one obtains 
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dQy 



(6) 



Q y v, 



xlm{[^ + e^(^)^; 



V Q0 



th 



-}• 



The last term inside the square bracket of this expres- 
sion gives no contribution when the electron trajectory 
is contained in a plane of specular symmetry of the crys- 
tal surface. 

As a test one can apply Eq. (6) to an electron moving 
in vacuum (e^ = ph = 1) at a distance zq from the surface 
of a non-magnetic medium described by e, in which case 



R, 



QO 



DPS 



P(W) 



0. Then, Eq. (6) reduces to 
2 f°° Tq dQ y 



-2r QZo 



x Im{ 



o 



r; 



Q +( ^ )2 



r d 



eT Q + T Q Tqc T q +T[ 



-}, 



Q 
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where T' Q = ^/Q 2 - ew 2 /c 2 , T Q = ^/Q 2 - lo 2 /c 2 , and 
the reflection coefficients of Eq. (6) have been taken 
from Fresnel's equations. 27 In the non-relativistic limit, 
= Tg = Q, so that the contribution of the reflec- 
tion of s waves vanishes and the loss probability becomes 
proportional to Im{— l/(e+ 1)}- 

III. ENERGY LOSS FOR SMALL LATTICE 
CONSTANTS 

When the crystal lattice constant is much smaller than 
the wavelength corresponding to a given frequency com- 
ponent u), the details of the crystal lattice cannot be re- 
solved by the electron, so that the medium behaves like 
a uniform material, characterized by an effective dielec- 
tric constant and magnetic permeability. For instance, a 
dilute system of spheres can be regarded as a set of inter- 
acting dipoles, which leads to the well-known Maxwell- 
Garnett formula. 28 For more compact systems the details 
of the microscopic structure becomes relevant via impor- 
tant multipolar interactions. 

These are the cases considered in Fig. 1 for a crystal 
composed of small aluminum spheres that are disposed 
in a simple cubic lattice of constant a = 5 nm. The 
solid curves of Fig. I (a) show the loss probability for an 
electron moving with velocity v — 0.4c parallel to the sur- 
face of a slab made up of six (100) layers of such crystal. 
The lost probability has been calculated by means of Eq. 
(6) for two different filling fractions / of the aluminum 
(6.5% and 30%, respectively). The dielectric function of 
aluminum has been approximated by a Drude expression 
e(uj) = 1 — lo 2 /lo(uj + ir]) with uj p = 15 cV and rj = 1 
eV. The wavelengths considered in the figure lie in the 
range A = 2-kc/u — 88.6 — 310 nm, so that we are in the 
A >> a limit. The maximum orbital quantum number 
used to achieve convergence on the figure is Z max = 6. 

Crystalline effects come from the interaction among 
spheres, and a method to determine their relative role 
consists in comparing these results with a calculation in 
which that interaction is suppressed. This is what the 
broken curves stand for. It is very clear that the inter- 
action among spheres is almost negligible for / = 6.5%, 
whereas it becomes very important at / = 30%. In the 
latter case, the sphere surfaces are closer to each other 
and their mutual electromagnetic coupling becomes rel- 
evant. The non-interacting case shows a peak at around 
uj\ = 8.7 cV that corresponds to the dipole Mie reso- 
nance in the small-sphere limit, given by the expression 
e(uji) — —2 (i.e., uj\ = w p /\/3)- The main effect of the 
sphere-sphere interaction (solid curve for / = 30%) con- 
sists in splitting this peak, in a similar way as splitting 
of degenerate levels (the Mie resonances here) occurs in 
atomic bonding. 

This is actually observed in the projected photonic 
band structure of this crystal, represented in Fig. 2 
for / = 30%. The structure is dominated by nearly 



fiat bands corresponding to localized excitations near 
the Mie resonances of the small isolated spheres, loi — 
w p \Jl/{l + 1), which lie in the 8.7-10.6 eV range (see the 
labels on the right hand side of the figure for the dipole 
and quadrupole Mie modes). The interaction between 
spheres gives rise to a complex structure that encom- 
passes regions of relative transparency outside that en- 
ergy range. The evanescent plane wave components of 
the external electron field are subject to the condition 
cj = Q v, which defines the straight lines shown in 
the figure. The loss spectrum can then be understood 
as originating from both absorption and direct coupling 
to propagating modes of the crystal. This last effect is 
clearly seen as a bump in the loss probability near 6.7 
eV, connected to the crossing of the noted straight lines 
with a low-energy region of allowed propagating modes. 

The dashed curves of Fig. 1(a) (uncoupled spheres) 
have been obtained in two different ways. The first one 
consists in calculating the reflection coefficients that ap- 
pear in Eq. (6) using the layer KKR method but ne- 
glecting the interaction among spheres within each layer 
and also the interaction among different layers. A sec- 
ond procedure consists in making use of the analytical 
expression for the energy loss probability of an electron 
near an isolated sphere (see Refs. [ 6] and [ 9, 10] for non- 
relativistic and fully-relativistic formulas, respectively) 
and summing over all impact parameters of the electron 
trajectory with respect to the spheres of the crystal. The 
results coming out of these two very different procedures 
cannot be distinguished on the scale of the figure, and 
this is a strong indication of the degree of convergence of 
our numerical calculations with respect to the number of 
plane waves used in the layer KKR method. 

These results do not support, however, previous cal- 
culations by Pendry and Martin-Moreno 19 for this exact 
system, where the number of features and their energy 
position for the filling fractions under consideration differ 
from ours. A possible lack of convergence of the trans- 
fer matrix method used by the authors for this three- 
dimensional system of metallic scatterers might be the 
reason of this discrepancy. In particular, some of their 
low-energy features around 6.5 eV could originate in the 
modes of the wedge associated to their space discretiza- 
tion procedure (resonances near that energy have also 
been found for the 90° aluminum wedge 9 ). 

The example offered in Fig. 1 illustrates what happens 
in the a << A limit, where it is reasonable to define 
an effective dielectric function e e s for the infinite crystal 
and to compare the results of the detailed, exact calcu- 
lation with those obtained for an electron moving par- 
allel to a homogeneous surface of a material described 
by such a dielectric function. Several recipes for defining 
e ff for granular materials have been given in the litera- 
ture, ranging from simple effective medium theories like 
Maxwell-Garnett's 28 to more elaborate ones that take 
into account the actual shape of the constituents, both 
for disordered composites 29 ' 30 ' 17 and for crystals. 29 ' 31-34 
Among the latter, one finds extensions of the Maxwell- 
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Garnett theory that go beyond dipolar interactions, ' 
spectral representations, 31 ' 33 or direct derivation of the 
light dispersion relation. 32 Here, we have used a differ- 
ent method that consists in finding the dielectric function 
e c ff of the equivalent homogeneous medium that has the 
same reflectance as the crystal surface; 26 the reflectance 
calculated for the crystal is well reproduced by Fresnel 
equations with a single value of e e s(u>) for each frequency 
u> within a 2% under the present conditions. The re- 
sults are represented by the solid curves of Fig. 1(b), 
where they are compared with Maxwell-Garnett's the- 
ory (broken curves). For / = 6.5% these two models 
are nearly identical, as Maxwell-Garnctt formula gives 
the dipole of the small isolated sphere correctly. How- 
ever, for / = 30% the deviation between both models 
is significant, since higher momenta are involved in the 
interaction between neighboring spheres . The magni- 
tude actually represented in the figure is the surface loss 
function, Im{— 1/(1 + e e ff)}, which is directly compara- 
ble to the loss probability of Fig. 1(a). The agreement 
in the position of the peaks between the detailed energy 
loss calculation and our effective medium theory is rea- 
sonably good, indicating that an effective homogeneous 
medium describes the solid appropriately within this en- 
ergy range at v = 0.4c; however, the relative weight of 
the features for / = 30% changes completely at lower 
velocities (dotted curve, v — 0.06c = 8.2 a.u.). 

Another example of a crystal that can be represented 
by an effective dielectric constant is given in Fig. 3, con- 
sisting of a simple cubic lattice with the same parameters 
as in Fig. 1, except that finite cylinders have been used 
instead of spheres. As a result, the medium is strongly 
anisotropic and characterized by different bulk plasmon 
modes when the electric field is directed parallel or per- 
pendicular to the cylinders [see Fig. 3(c)-(d)]. The loss 
probability [Fig. 3(a)] is shown to share most of the fea- 
tures of the surface loss function for anisotropic media 
[Fig. 3(b)], except for the 7.8 eV peak in the loss spec- 
trum with / = 30%, which might be connected to the 
proximity effect when the electron starts sensing the non- 
uniform character of the surface via evanescent waves. 



IV. ENERGY LOSS FOR LATTICE CONSTANTS 
COMPARABLE TO THE WAVELENGTH 

For lattice constants comparable to the wavelength as- 
sociated to the energy losses under consideration, one 
can no longer define the effective dielectric constant of 
an equivalent homogeneous medium. Then, it is useful 
to relate the loss spectra directly to the photonic band 
structure. This has been done in Fig. 4 for a crystal of 
aluminum spheres immersed in a dielectric with e = 3. 
The features of the loss spectra [Fig. 4(c)] are strongly 
correlated with the band structure [Fig. 4(a)-(b)]. In par- 
ticular, the pseudo-gaps near 3.2 eV are translated into 
a dip in the loss probability. 



Using electrons to analyze the crystal permits explor- 
ing regions of the band structure that are not accessi- 
ble to external light, but that can be reached via the 
evanescent waves contained in the perturbing field of the 
electrons. This is the case of Fig. 4 for v/c = 0.5 [solid 
curve in Fig. 4(c)], well below the Cherenkov threshold 
v/c= 1/ V3, that permits to observe how the peak in the 
loss spectrum near 2.7 eV results from the coupling with 
propagating modes of the crystal as a result of the in- 
tersection of the external field components [straight solid 
lines in Fig. 4(a)] and the allowed regions of propagation 
(shaded areas). 

At larger velocities [v/c = 0.75, dashed curve in Fig. 
4(c)], Cherenkov radiation would be produced in the ab- 
sence of the crystal, that would result in a loss probabil- 
ity independent of uj (see little arrow near the horizontal 
axis of the figure). This is clearly seen in the big over- 
lap of the external field with the regions of propagation 
within the crystal inside the low-energy region results in 
an enhanced loss probability as compared to v/c = 0.5. 
However, the presence of the crystal modulates the loss 
spectrum. 

Part of the energy lost by the electron must be con- 
verted into the so-called Smith-Purcell radiation. 21,10 In 
particular, when transparent materials are used to build 
the crystals, the light emission probability must coincide 
with the energy-loss probability. For a crystal of finite 
thickness, like the one considered in Fig. 5, consisting of 
8 (lll)-layers of an inverted Si opal (e = 11.9) with a fill- 
ing fraction of 67%, part of this emission occurs towards 
the side of the crystal opposite to the electron trajectory 
(here, the electron is taken to be moving parallel to the 
crystal surface). The intensity of the emitted light [Fig. 
5(c)] has been calculated from the integral of the Poynt- 
ing vector far from the crystal (see Ref. [ 23] for more 
details), and it presents a strong dip near a/A w 0.75, 
where a = 1220 nm is the lattice constant. The region 
of emission depletion is actually contained within a full 
band-gap of the infinite crystal [see Fig. 5(a)]; this is also 
seen in the transmission of light incident on the crystal 
both normal to the surface [dashed curve in Fig. 5(b)] 
or with an angle corresponding to the Cherenkov radia- 
tion produced by the electron in the plane defined by the 
trajectory and the surface normal (solid curve). 

V. ENERGY LOSS IN THE BULK OF A 
CRYSTAL 

So far we have considered an electron moving in front 
of the surface of a photonic crystal. When the electron is 
moving in the bulk of an infinite crystal, the electric field 
can also be written in terms of reflectance coefficients if 
the trajectory is contained in a low-index plane that does 
not intersect any of the crystal objects. The reflectance 
in question is that of the lower and upper semi-infinite 
crystals into which the entire crystal is divided by the 
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noted plane. The corresponding reflectance matrices will 
be denoted R\ and R2, respectively. They contain the 
amplitudes of reflected beams, labeled by vectors G of 
the 2D reciprocal lattice of the plane under consideration. 
In particular, for a fixed value of Q, i?i GG / is given by 

^q+G'.g-G'i as defined by Eq. (1). Similarly, i?f GG , 
connects Kq +G wave components with waves reflected 
from the upper surface with momentum Kq +G , . 

This is represented schematically in Fig. 6, where the 
reflectance matrices of the lower and upper semi-infinite 
crystals are R\ and R2, respectively, and upwards (down- 
wards) arrows represent waves of momentum Kq +g 

(Kq +g ). The upper and lower surfaces have been sepa- 
rated in the figure for the sake of clarity, but they will be 
considered to lie on the same plane (a plane that contains 
the electron trajectory) in what follows. 

Using Eq. (5) for the loss probability and a straightfor- 
ward extension of Eq. (4) for the reflected electric field 
including all terms of Fig. 6, the loss probability per unit 
of path length averaged over all impact parameters par- 
allel to the crystal surfaces is found to be 



P r (w) = —J 
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where the matrices ^4 G CT G /, B G a G ,, C G a G ,, and D G a G , are 
defined in Fig. 6, and only the elements G = G' = 
enters this expression. 

The above formula has been applied to calculate the 
energy loss probability for an electron moving in between 
two parallel crystal slabs of 32 layers each surrounded by 
Si, as shown in Fig. 7. The crystal parameters are the 
same as in Fig. 5(c), so that an absolute band gap shows 
up near a/ A = 0.75. The electron energy (100 keV) is 
large enough to produce Cherenkov light, which must 
be confined in between the two crystals for wavelengths 
lying in the gap region. This is actually observed as a 
pronounced dip, connected to the fact that no electro- 
magnetic modes can be created that escape through the 
crystals, except for a small transmission due to their finite 
thickness. However, some radiation can escape though 
modes that are trapped in the slab formed by the two 
crystals, so that the loss probability does not actually 
reach a zero value in the dip. 



VI. CONCLUDING REMARKS 

The energy loss of fast electrons interacting with pho- 
tonic crystals has been calculated for dielectric, metallic, 
and metalo-dielectric systems. 

When the wavelength of the radiation associated to the 
energy loss is much larger than the lattice constant, the 
crystal can be regarded as a continuous medium, charac- 
terized by an effective dielectric function, as it has been 
shown for crystals of small aluminum spheres in Fig. 1 
and small finite cylinders in Fig. 3. The effective medium 
is highly anisotropic in the latter case. The effective di- 
electric function has been calculated in both cases from 
the reflectance coefficient of the crystal, and this has been 
shown to contain most of the information needed to un- 
derstand the calculated energy-loss spectra. 

When the lattice constant is comparable to the wave- 
length, the idea of an equivalent effective continuous 
medium is no longer valid, and one has to rely on the 
detailed band structure of the crystal to understand the 
loss spectra [Figs. 4, 5, and 7]. 

Finally, the interaction of the electron with the crys- 
tal produces Smith-Purcell radiation, which contributes 
to the energy loss. If the crystal is composed of non- 
absorbing materials, the energy loss probability and the 
light emission probability must coincide, as shown in Fig. 
5 for an inverted Si opal, in which case the light emitted 
after transmission through the crystal shows dips that 
are directly connected to the presence of photonic band 
gaps. For electrons moving in a region surrounded by 
photonic crystals a dip is also observed in the loss prob- 
ability within the gap energies (Fig. 7). 

Our hope is that the present work can provide a stim- 
ulus to use fast electrons in the analysis of photonic crys- 
tals as a way to bring a source of evanescent radiation 
(the electromagnetic field of the electron in vacuum) into 
close contact with the crystal and also to probe regions 
of the crystal that would not be easily accessible to other 
sources of external electromagnetic radiation. 
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FIG. 1. (a) Energy loss spectra for an electron moving par- 
allel to the [100] direction of the (100) surface of a simple cubic 
crystal made up of six layers of aluminum spheres in vacuum 
with lattice constant 5 nm and two different filling fractions 
(see labels). The electron is moving at a distance of 1 nm from 
the sphere surfaces with a velocity v = 0.4c (solid and dashed 
curves) and v — 0.06c (dotted curve). The dashed curves 
are obtained by neglecting the interaction among spheres, 
whereas the solid curves and the dotted curve correspond to 
the full solution of Maxwell's equations, (b) Surface loss func- 
tion for the same crystal as in (a) using the effective dielec- 
tric function obtained from the Maxwell-Garnett expression 
(broken curves) and from the multiple-scattering method ex- 
plained in the text (solid curves). 
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FIG. 2. Photonic band structure of one of the crystals con- 
sidered in Fig. 1, consisting of a simple cubic lattice of alu- 
minum spheres in vacuum with a lattice constant of 5 nm 
and a filling fraction of 30%. The figure shows the band 
structure projected on the (100) surface. The horizontal axis 
represents the excursion along the points indicated in the 
inset within the parallel momentum plane. The expression 
e(a>) = 1 — uJp/cj 2 has been used for the aluminum dielec- 
tric constant with u> p = 15 eV. When the damping is taken 
as ij — > + , one obtains regions of allowed electromagnetic 
propagation, that is, combinations of the energy and the par- 
allel momentum components for which some eigenstates have 
a real perpendicular momentum. They define the shaded re- 
gion in (a). However, for a realistic value of the damping 
parameter (77 = 1 eV), all perpendicular momentum compo- 
nents are complex and they represent evanescent waves within 
the crystal. The contour plot in figure (b) represents the low- 
est value of the imaginary part of the perpendicular momen- 
tum multiplied by the lattice constant (i.e., the minimum of 
Im{fc z a}) for all eigenstates with a given parallel momentum. 
The boundaries of the shaded areas of (a) are shown in (b) as 
continuous curves. 
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FIG. 3. (a) Energy loss spectra for an electron moving par- 
allel to the [100] direction of the (100) surface of a simple cubic 
crystal made up of six layers of aluminum cylinders in vacuum 
with lattice constant 5 nm and two different filling fractions 
(see labels) . The height of the cylinders is equal to the diam- 
eter in all cases. The electron is moving at a distance of 1 nm 
from the cylinder surfaces with a velocity v = 0.4c. (b) Sur- 
face loss function for the same crystal as in (a) using the effec- 
tive dielectric function obtained from the multiple-scattering 
method explained in the text, (c) Anisotropic dielectric func- 
tion for the crystal considered in (a) with a filling fraction of 
the aluminum of 6.5%. (d) Same as (c), for a filling fraction 
of 30%. 
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FIG. 4. (a)-(b) Photonic band structure projected on the 
(100) surface and represented as in Fig. 2(a)-(b) for an fee 
lattice of aluminum spheres surrounded by a medium of di- 
electric function equal to 3. The filling fraction of the spheres 
is 20% and the lattice constant is a = 105.2 nm, comparable to 
the wavelength A (see scale on the left hand side of the figure) . 
(c) Electron energy loss spectra for an electron moving with 
velocity v = 0.5c (solid curve) and v = 0.75c (broken curve) 
at a distance of 3.2 nm from the surface of the outer-most 
spheres of a semi-infinite crystal like that considered in (b). 
The trajectory is directed along the [100] direction (TX). The 
frequency and parallel momentum of the plane wave compo- 
nents of the external electron field are related by the condition 
ui — Q ■ v, which defines the straight lines shown in (a) and 
(b). 




FIG. 5. (a) Photonic band structure projected on the (111) 
surface for an fee lattice of air spherical voids in Si (e = 11.9). 
The radius of the spheres is 0.342 times the lattice constant, 
a. (b) Dashed curve: transmittance of light incident nor- 
mal to 8 layers of the crystal considered in (a). Solid curve: 
transmittance for an angle of incidence corresponding to the 
Cherenkov angle 6 C = 68.1° for a 100-keV electron moving 
parallel to the same crystal. The crystal is surrounded by Si 
on both surfaces . (c) Probability of emitting light on the 
opposite side of the crystal with respect to the electron tra- 
jectory. The distance from the latter to the sphere surfaces is 
876 nm and the lattice constant is a = 1220 nm. The proba- 
bility is given per unit of path length in microns and per unit 
of emitted-photon wavelength also in microns. The region of 
frequency and momentum where the field of the electron in- 
side a homogeneous infinite Si medium takes non-zero values 
is represented by thick solid lines in (a) . 
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FIG. 6. Schematic representation of the construction of the 
electric field in the bulk of a photonic crystal in terms of the 
reflectance of two semi-infinite crystals and definition of the 
matrices A, B, C, and D used in Eq. (7). See Sec. V for more 
details. 
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FIG. 7. Energy loss probability for a 100-keV electron trav- 
eling parallel to two photonic crystals with the same param- 
eters as in Fig. 5(c). The spacing between the two crystals 
is 300 nm, which is filled with Si (this is exaggerated in the 
inset for the sake of clarity). The electron is traveling at the 
same distance from the two crystals. 
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